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ABSTRACT 

Exploiting the relationship between CO luminosity and line width determined for fainter starbursts, we identify 
and remove lensed sources from the widest Herschel extragalactic imaging survey to yield a sample of intrinsi- 
cally luminous galaxies. Panchromatic observations of the best candidate HyLIRG then led to the discovery of 
at least four luminous galaxies across a ^100-kpc region at z = 2.41. Our technique separates HyLIRGs from 
the more numerous, strongly lensed population, and can also pinpoint clusters of starbursting proto-ellipticals. 
Via sub-arcsecond interferometric imaging we have measured accurate gas and star-formation surface densi- 
ties. The two brightest galaxies span ~3 kpc FWHM in submm/radio continuum and CO /=4-3, and double 
that in CO / = 1-0. The broad CO line is due partly to the multitude of constituent galaxies and partly to 
large rotational velocities in two counter-rotating gas disks - a scenario predicted to lead to the most intense 
starbursts, which will therefore come in pairs. The disks have Mdyn of several x 10^^ M©, and gas fractions of 
~ 40%. Velocity dispersions are modest so the disks are unstable, potentially on scales commensurate with 
their radii: these galaxies are undergoing extreme bursts of star formation, not confined to their nuclei, at close 
to the Eddington limit. Their specific star-formation rates place them > 5 x above the main sequence, which 
supposedly comprises large gas disks like these. Their high star-formation efficiencies are difficult to reconcile 
with a simple volumetric star-formation law. N-ho&y and dark matter simulations suggest this system is the 
progenitor of a B(inary)-type ^ lO^^-^-M© cluster. 

Keywords: galaxies: high-redshift — galaxies: starburst — submillimeter: galaxies — infrared: galaxies — 
radio continuum: galaxies — radio lines: galaxies 
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L INTRODUCTION 

Of the known denizons of the galaxy zoo, hyperlumi- 
nous infrared (IR) galaxies (HyLIRGs, Ljr > 10^^ L©, where 
the IR luminosity is measured across Arest = 8-1000 /im) are 
amongst the rarest and most extreme. They provide excellent 
laboratories with which to confront the most recent hydro- 
dynamic simulations of isolated and merging galaxies (e.g. 
iHayward et al. 201 1). The Ljr of a HyLIRG implies a stagger- 
ing st ar-formation rat e, SFR > 1000 M© yr"^ (Kennicuttll998L 
for a iChabrieri 120031 initial mass function, IMF), unless the 
IMF is top-heavy, or there is a substantial contribution to Ljr 
from a deeply obscured AGN, in which case we are seeing 
intense star formation accompanied by the rapid growth of a 
massive black hole ( Alexander et al. 2008). Either way, we 
are witnessing galaxy formation at its most extreme. 

How best can we identify the most luminous star-forming 
HyLIRGs? A promising method involves searching amongst 
the brightest sources detected in the widest surveys with Her- 
schel. These comprise blazars, low-redshift spirals and some 
intri nsically fainter sour ces that have been strongly lensed 
(e.g. ' Negrello et al.l 1201 0). However, there is also the in- 
triguing possibility that such samples may contain unlensed 
and thus intrinsically luminous galaxies - the rarest HyLIRGs 
are predicted to have a space density, ^ 10~^Mpc~^ and 
may mark the sites of today's clusters (|Negrello et al.l 120051: 
iLapi et aLllMTli) . 

Of the facilities providing redshifts for bright Herschel 
sources, the 100-m Robert C. Byrd Green Bank Telescope 
(GET), together with the ultrawide-bandwidth Zpectrome- 
ter cross-correlation spectrometer and Ka-ban d receiver, have 
been amongst the most effective (e.g. Swinbank et al.ll20TQl: 
iFrayer et al]l201 ll) . Having added consi derably to the sam- 
ple of such sources with known redshifts, iHarris et aP (|2012|) 
identified a flat trend of L^q with ^^C^^O J= 1-0 linS width. 
This contrasted with the steep pow er-law rel ation between 
L[^Q and CO line width found by Both well et al. (2013) for un- 
lensed SMGs, as expected if the most gas-rich gal axies tend to 
live in the most massive gravitational potentials. iHarris et aH 
argued that intrinsically fainter sources require a higher tens- 
ing magnification, /i, to rise above the observational detection 
threshold for CO / = 1-0, which is approximately constant 
in flux. The flat trend seen in the GET sample is then best 
interpreted as tensing magnification acting on a population 
with intrinsically steep number counts. Those GET sources 
closest to the power-law fit seen for SMGs - typically those 
with the widest lines - are likely to be suffering the least tens- 
ing magnification. This is where we might expect to find any 
HyLIRGs that are lurking amongst the lensed starbursts. 

Our approach here, therefore, is to search for HyLIRGs 
amongst those bright Herschel lens candidates with the 
broadest CO lines, starting with th e best exa mple in 
the sample observed with GET by ^ Harris et al.i. HAT- 
LAS J084933.4+021443 (hereafter HATLAS J084933, with 
350-/im flux density, ^'350 = 293 mJy), which displays a line 
consistent with CO / = 1-0 at zlsr = 2.410 ±0.003, with a 
fuU width at half maximum (fwhm) of 1180 ± 320 km s"^ 
and ^coi-o = 0.83 ± 0. 19 mJy. The amplification predicted by 
IHarris et al] for HATLAS J084933 is consistent with unity, 
/i = 2±l. 

In the next section we describe an extensive set of ob- 
servations. We present, analyse, interpret and discuss our 

Hereafter, CO refers to ^^C^^O unless stated otherwise. 



reduced images, spectra and cubes in §|3] finishing with 
our conclusions in ^ We adopt a cosmology with Hq = 
71 kms"^ Mpc"\ Qrn = 0.27 and Qa = 0.73, so I" equates to 
8.25 kpc at z = 2.41. 

2. OESERVATIONS AND DATA REDUCTION 

Herschel imaging, unde rtaken as part of t he wide-field H- 
ATLAS imaging survey (lEales et al.l 120101) . led to the se- 
lection of HATLAS J084933 as a potentially distant, lensed 
starburst: distant, because it is an example of a so-called 
'350-/im peaker', with its thermal dust peak between the 
250- and 500-/im bands; lensed, because models of the far- 
IR/submm source counts suggest that the majority of objects 
with 5*500 > lOOmJy are expected to be either local (z < 0.1) 
or lensed (e.g. N egrello et al.ll2010|) . along with the occasional 
flat- spectrum radio quasar. 

On the basis that it represented the best chance of finding an 
intrinsically luminous system rather than a lensed galaxy, fol- 
lowing the arguments laid out in HATLAS J084933 was 
selected for observations with the Jansky Very Large Array 
(JVLA). Here we detail the JVLA observations and those ob- 
tained with other facilities that followed as a consequence of 
our initial findings, presenting these datasets in order of de- 
creasing wavelength. 

2.1. JVLA CO J =1-0 and 5 -GHz continuum imaging 

Whilst the National Radio Astronomy Observatory's 
(NRAO's) JVL43 was in its DnC, C, E, EnA and A configu- 
rations, between 2012 January and 2013 January, we acquired 
^30 hr of Ka-band data, scheduled dynamically to ensure ex- 
cellent atmospheric phase and pointing stability. We recorded 
two sets of eight contiguous baseband pairs, 1024 x 2-MHz 
dual-polarization channels in total. We tuned the first set of 
basebands to cover ^^C^^O and ^^c^^O / =l-0. The ^^C^^O 
/= 1-0 transition (i^rest = 115.271203 GHz.l Morton & Noreaul 
IT994I) was placed in the second set of baseband pairs, offset- 
ing down by 64 MHz to 33.740 GHz to avoid baseband edges. 

Around 4 hr of pseudo-continuum A-configuration C-band 
data were also obtained, recording 1024 x 2-MHz dual- 
polarization channels every 1 sec across 4.2-6.5 GHz, with 
a small gap. Typically, 0.3 GHz was lost to severe radio- 
frequency interference near 6.1 GHz, yielding a band center 
of 5.1GHz (5.9 cm). 

PKS J0825+0309 and PKS J0839+0319 were observed ev- 
ery few minutes to determine accurate complex gain solutions 
and bandpass corrections in the Ka and C bands, respectively. 
3C 286 was observed to set the absolute flux density scale, 
and the pointing accuracy was checked locally every hour. 

The data were reduced using AIVS. The basebands were 
knitted together using the NOIFS and VBGLU tasks, yielding 
uv datasets with two intermediate frequencies, each compris- 
ing 512 X 2-MHz channels. 

For the spectral-line data, these were then imaged in groups 
of four channels (71 kms~^), with natural weighting. A vari- 
ety of Gaussian tapers were used to weight the data, with dis- 
tances to the 30 per cent point of the Gaussian ranging from 
100 to 1,000 kA, to form a number of cubes with FWHM spa- 
tial resolutions ranging from 0.53'' x 0.51'' to 2.0" x 1.8". 
The line-integrated map shown in Fig. [T]has a 1.1" x 1.0" 
synthesized beam. 

This work is based on observations carried out with the JVLA. The 
NRAO is a facility of the NSF operated under cooperative agreement by As- 
sociated Universities, Inc. 
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The naturally- weighted 5.1 -GHz pseudo-continuum data 
were imaged, averaging al l of the channels, loosely f ollowing 
the techniques outlined by lQwen & MorrisonI (I2QQ8I) . The re- 
sulting map has a spatial resolution of 0.53'' x 0.42'' (PA, 26°) 
and an r.m.s. noise level of 2.9 /iJy beam~^ (Fig.[T]). 

2.2. CO J =3-2 imaging from CARMA 

We used CARM43 to observe the CO 3-2 line (z/rest = 
345.795991GHz, redshifted to z/obs= 10 1.406 GHz). Obser- 
vations were carried out with the 3 -mm receivers and the 
CARMA spectral-line correlator, with an effective band- 
width of 3.7 GHz per sideband and a spectral resolution of 
15kms~^ Four tracks were obtained between 2012 Febru- 
ary 13-17 in the C configuration ( 15-3 5 2-m baselines), with 
14 usable antennas and a total on-source observing time of 
13.5 hr. The nearby quasar, PKS J0757+0956, was used for 
complex gain calibration. The bandpass shape and abso- 
lute flux calibration (the latter good to ~15%) were derived 
from observations of Mars, PKS J0854+2006 and 3C 84. Us- 
ing natural weighting, we obtained a FWHM beam size of 
2.1" X 1.8" (PA, 126°), with a noise level, cr = 0.3 mJybeam-i 
(Fig.D. 

2.3. CO J =4-3 imaging from IRAMPdBI 

During 2012 February we obtained data using six 15-m an- 
tennas in the most extended configuration of the Institut de 
Radioastronomie Millimetrique's Plateau de Bure Interferom- 
eter (IRAlvQPdBI), with baselines ranging up to 800 m. 

The observing frequency was set to 135.203 GHz, the 
redshifted frequency of the CO / = 4 - 3 line (i^rest = 
461.04077 GHz). The weather conditions were excep- 
tionally good with a precipitable water vapour of around 
2 mm and 7;ys - 100 K on average. PKS J0825+0309 
and PKS J0909+0121 were used to calibrate the complex 
gains, while 3C84 and MWC349 were chosen as band- 
pass and absolute flux calibrators. The r.m.s. noise level 
is 0.97mJybeam~^ in 4-MHz- wide spectral channels, for a 
synthesized beam measuring 1.0" x 0.5" FWHM (PA, 24°). 
Using only the line-free channels, the r.m.s. noise level is 
47 /iJy beam~^. 

We also report continuum photometry at 1 .36 mm from ob- 
servations obtained throughout 2012 as part of a search for 
water that will be reported elsewhere. 

2.4. SVO-fim continuum imaging from the SMA 

Approximately 1.4, 4.9 and 1.7 hr of integration were ob- 
tained using the Submillimeter Array (SMAH) in its com- 
pact, extended and very-extended array configurations during 
2011 December 8, 2012 January 31 and 2012 March 31, re- 
spectively (2011B-S044). The receivers were tuned such that 
the upper sideband was centered at 345 GHz (870 /im), with 
8 GHz of single-polarization bandwidth in total. 3C 84 was 

CARMA construction was derived from the states of Maryland, Califor- 
nia and Illinois, the James S. McDonnell Foundation, the Gordon and Betty 
Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the Uni- 
versity of Chicago, the Associates of the California Institute of Technology, 
and the NSF. Ongoing CARMA development and operations are supported 
by the NSF under a cooperative agreement, and by the CARMA partner uni- 
versities. 

41 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 

4^ The SMA is a joint project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy and Astrophysics 
and is funded by the Smithsonian Institution and the Academia Sinica. 



used as the bandpass calibrator and Titan was used for abso- 
lute flux calibration. PKS J0825+0309 and PKS J0909+0121 
were again used to track and check the phase and gain. 

A ROBUST = 2 weighting scheme resulted in a 0.83'' x 
0.67'' beam FWHM (PA, 86°) and a noise level, a = 
0.77mJybeam-i (Fig.Q). 

2.5. 350-iim continuum imaging with CSO/SHARC-2 

A total of 80min of data were collected at 350 /im dur- 
ing periods of excellent weather between 2012 January 12- 
13 at the 10.4-m Caltech Sub millimeter Observat ory (CSO), 
using the SHARC-2 camera (iDowefl et al.ll2003h . We took 
advantage of the open-loop, actuated surface of the telescope 
(Leong et al. 2006) and we checked and updated the focus 
settings several times each night, resulting in near-Gaussian, 
diffraction-limited (8.5 " FWHM) be am profiles. Data were 
reduced using CRUSH (lKovacsll2008|) and the resulting maps 
have a noise level, cr = 8 mJy beam~^ . 

2.6. Continuum imaging and FTS spectroscopy from 

Herschel 

The acquisition and reduction of Her schel (iPilbratt et al.l 
'20T0^ parallel-mode data from SPIRE (Griffi n etal.lloTl ^ 
and PACS ( Poglitsch et al. 2010) for the 9-hr Science Demon- 
stration Phase field of ATL AS (Eales et al. 2010) are de- 
scri bed in detail by llbar et all (l2Q10a) . Pascale et al.l (1201 Ih 
and iRigbv et al.l (1201 Ih . Here, we use the original SPIRE 
imaging and have obtained much deeper PACS data from the 
OTl programme, 0T1_RIVIS0N_1. We recorded data simul- 
taneously at 100 and 160 /im whilst tracking across the target 
at 20" s~^ in ten 3-arcmin strips, each offset orthogonally by 
4". A total of 180 s on source was acquired for each of two 
near-orthogonal scans. These data were tackled with a vari- 
ant of the pipeline developed by Ibar et al. ( 2010b), reaching 
G ^ A and 7mJy at 100 and 160 /im, respectively. The 160- 
/im PACS image is shown in Fig.[T] 

Spectroscopic data covering Aobs = 1 97-670 /im were ac- 
quired during 2012 November using the SPIRE Fourier Trans- 
form Spectrometer at the highest available spectral resolution 
as part of ot2_rivison_2. These data comprise 100 forward 
and 100 reverse scans, totalling almost 14ksec, with the cen- 
tral pixel centered on the SPIRE centroid, as reflected by the 
target's full lAU nomenclature. 

2.7. Infrared continuum imaging from Spitzer, VISTA and 
the Hubble Space Telescope 

3.6- and 4.5-/im imag es were acquired using the Infrared 
Array Camera (IRAC - Fazio et al.l l2004h aboard SpitzeJ^ 
(Werner etal. 2004) on 2012 June 24 as part of program 
80156. The imaging involved a 36-position dither pattern, 
with a total exposure time of just over 1 ks, reaching r.m.s. 
depths of 0.3 and 0.4 /iJy at 3.6 and 4.5 /im, respectively. Cor- 
rected basic calibrated data, pre-processed by the Spitzer Sci- 
ence Center, were spatially aligned and combined into mo- 
saics (shown in Fig. [T]) with a re-sampled pixel size of 0.6" 
and ang ular resolution of 2-2.5", using version 18.5.0 of 
MOPEX (iMakovoz & Marleaull2005h . 

Images were obtained in Z,7,/,^ and as part of 
VIKING, a public survey with the 4-m Visible and Infrared 

This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, CaHfornia 
Institute of Technology under a contract with NASA. 
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Survey Telescope for Astrono my (VISTA) at Paranal, Chile 
(lEmerson & Sutherlandll2QlQk . with an image quality of 0.9 
FWHM, as shown in Fig.[Tl 

A SNAPshot observation was obtained with the Hub- 
ble Space Tele scop (HST) on 2012 February 08, as part 
of Cycle- 19 proposal 12488, using Wide-Field Camera 3 
(WFC3) with its wide / filter, FllOW. The total exposure 
time was 252 s. Data were reduced using the IRAF MultiDriz- 
zle package. Individual frames were corrected for distortion, 
cleaned by cosmic rays and other artifacts and median com- 
bined. The resulting ^2' x2' image was re-sampled to a finer 
pixel scale of 0.064'' and is shown in Fig.[T] 

2.8. Optical spectroscopy from Keck 

We observed several optical/IR sources in the vicinity of 
HATLASJ084933 using the DEIMOS spectrograph on the 
10-m Keck-ii telescopS during 2012 March 02. DEIMOS 
was used with its 600ZD grating (ruled with 600 lines mm~^), 
the GG455 filter and a 1'' slit at a central wavelength of 
720 nm. Two positions were targeted, each with a different 
PA, for 40 and 75 min. 

3. RESULTS, ANALYSIS AND DISCUSSION 

3.1. Basic morphological description 

At the spatial resolution of SPIRE, HATLASJ08493 is 
essentially unresolved. However, with only modestly im- 
proved spatial resolution, via deep imaging with PACS at 
100-160 /im and CSO/SHARC-2 at 350 /im, the system is 
resolved into two sources, referred to hereafter as W and T 
(Fig.ID. 

Moving to higher spatial resolution, via interferometric 
imaging with SMA, IRAM PdBI and CARMA at rest-frame 
255-880 /im, we see a consistent morphological picture for 
the emission from cool dust in the HATLAS J084933 sys- 
tem (Fig. [T]). The dust continuum emission is dominated 
by W and T, which are separated by 10.5'', or ~85kpc in 
the plane of the sky, with W the brightest of the two from 
the optical all the way through to the radio regime. W is 
amongst the brightest unlensed SMGs ever observed inter- 
ferometrically at submm wavelengths, with T not far behind 
fcear et al. 2000; Lutz et al. 2001; Younser et al. 2007, 2008, 
[20091: [Wang et al.l 120071: lAravena et al.l I20ia ISmolcic et al.l 
l2012allb|) . Amongst the unlensed SMG population, its flux 
density has been equalled only by SMM J12371 1+622212 
(also known as GN20, = 22.9 ±2.8 - iPope et al.1 120051: 
llono et al.l[2QQ6h . 

Aided by our interferometric astrometry we can see that at 
optical/IR wavelengths the westernmost clump, T, is barely 
visible as a red extension to a lenticular galaxy (Fig. [TJ. W 
coincides with the extremely red component of a red/blue pair, 
separated by ^2". Based on experience with SMGs such as 
SMM J02399-0136 (e.g. Ferkinhoff et al, in preparation), we 
expected to find evidence of an AGN towards the compact 
blue source, but found instead that it is a star - our optical 

Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with pro- 
gram 12488. 

Data presented herein were obtained at the W. M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute of 
Technology, the University of California, and NASA. The Observatory was 
made possible by the generous financial support of the W. M. Keck Founda- 
tion. 



Table 1 

Continuum flux densities 







W 


T 




0.88 


/im 


8.16±1.37 


5.80 ±0.98 


/iJy; VISTA Z 


1.02 


/im 


10.5 ±1.8 


8.65 ±1.45 


/iJy; VISTA Y 


1.1 


/im 


13.6±2.3 


11.2±1.9 


/iJy; FllOW 


1.25 


/im 


16.8±2.8 


13.0±2.2 


/iJy; VISTA / 


1.65 


/im 


23.5 ±4.0 


19.2±3.2 


/iJy; VISTA// 


2.15 


/im 


38.1 ±6.4 


26.1 ±4.4 


/iJy; VISTA 


3.6 


/im 


67.1 ±18.6 


47.7 ±17.4 


/iJy; IRAC 


4.5 


/im 


76.6 ±20.1 


45.8 ±15.2 


/iJy; IRAC 


100 


/im 


23 ±7^ 


19 ±7^ 


mJy; PACS 




/im 


91 ±12^ 


40 ±12^ 


mjy, rAUo 


250 


/im 


242 ±18^ 


mJy; SPIRE 


350 


/im 


293 ±22^ 


mJy; SPIRE 


500 


/im 


231 ±19^ 


mJy; SPIRE 


870 


/im 


25±2 


19±2 


mJy; SMA 


1.36 


mm 


8.3 ±0.5 


7.5 ±0.5 


mJy; PdBI 


2.21 


mm 


1.1±0.1 


0.8±0.1 


mJy; PdBI 


3.0 


mm 


0.33±0.14 


3a < 0.42 


mJy; CARMA 


8.8 


mm 


0.043 ±0.010 


3cr < 0.030 


mJy; JVLA 


5.9 


cm 


0.177±0.015 


0.065 ±0.015 


mJy; JVLA 






C 


M 




0.88 




5.96 ±1.00 




//Tv VISTA 7 


1.02 


/im 


5.46 ±0.92 




/iJy; VISTA Y 


1.1 


/im 


7.66 ±1.29 




/iJy; FllOW 


1.25 


/im 


5.72 ±0.96 




/iJy; VISTA/ 


1.65 


/im 


5.09 ±0.86 




/iJy; VISTA// 


2.15 


/im 


11.1±1.9 




/iJy; VISTA /Ts 


3.6 


/im 


19.5±5.9 




/iJy; IRAC 


4.5 


/im 


17.7±5.3 




/iJy; IRAC 


870 


/im 


4.6±1.9 


6.9±1.6 


mJy; SMA 


1.36 


mm 


1.8±0.2 


1.5±0.2 


mJy; PdBI 


2.21 


mm 


0.25 ±0.06 


0.31 ±0.08 


mJy; PdBI 


5.9 


cm 


0.022 ±0.006 


0.020 ±0.005 


mJy; JVLA 



^ Errors determined by placing many apertures across the map, with 3- and 
5-% calibration uncertainties added in quadrature at 100 and 160 /im, respec- 
tively. 

^ Errors include the contribution due to confusion and a 7-% calibration un- 
certainty has been added in quadrature (Valiante et al., in preparation). 



spectroscopy reveals absorption due to Mg as well as the Na 
D doublet. 

There is no obvious trend with wavelength for the relative 
IR/submm flux density contributions from W and T (listed 
in Table [B, Sw/St = 1.22 ±0.33, 2.22 ±0.42, 1.27 ±0.19, 
1.32 ±0.17, 1.13 ±0.10 and 1.23 ±0.25 at 100 /im, 160 //m, 
350 /im, 870 /im, 1.36 mm and 2 mm, respectively. The 
error- weighted mean for the IR/submm wavelength regime is 
1 .22 ± 0.07. Moving long ward to rest- frame 1 .7 cm, where we 
expect the emission to be dominated by synchrotron radiation 
related to supernova remnants (SNRs), HATLAS J084933 is 
also dominated by component W (Fig. [T]), with S^/Sj = 
2.7±0.7. 

Continuum emission from dust is also seen east of W, from 
a galaxy we label M in Fig.[T] M is detected securely in 870- 
/im, 1.36-mm, 2.21-mm and 5.9-cm continuum (Table[T]). 

Still further to the east, emission can be seen from a very 
red, morphologically distorted galaxy - labelled C in Fig. [l] 
- which bears a strong resemblence to the Antennae (e.g. 
IWhitmore & Schweized [T99l iKlaas et al.l r2010h . C is de- 
tected weakly in 870-/im, 1.36-mm, 2.21-mm and 5.9-cm 
continuum (Table [T]). 

As we shall see in more detail later ( ^3.5> , C and M are also 
detected in our CO / = 1-0 and /=4-3 spectral-line imaging, 
i.e. they also lie at z = 2.41, alongside W and T. 

Both M and C are coincident with IRAC 3.6-4.5-/im emis- 
sion, suggesting that these are not tidal tails or gaseous 
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Figure 1. Morphology of HATLAS J084933 at moderate-to-high spatial resolution. Top left: 160- and 870-/im continuum emission from dust, as detected by 
PACS and the SMA, superimposed on a three-color representation of the VISTA z,/ and data. We see two bright components: W is centered on a red galaxy, 
near a blue star; T lies at the end of a lenticular galaxy at z = 0.3478, coincident with its red southern tip. Dust emission can also be seen coincident with two 
further galaxies, labelled C and M, to the east of W. Both C and M can also seen in continuum at 5.9 cm (see lower right), at 1.36 mm, 2.21 mm, and in CO 
/= 1-0 and / = 4-3. Top right: CO / = 3-2 emission superimposed on a three-color representation of the Spitzer IRAC 3.6- and 4.5-/im imaging and a heavily 
smoothed VISTA J+H+K^ image as the blue channel. CO emission is coincident with the two submm-bright clumps, W and T. An apparent extension south of T 
is mirrored in the IRAC image, suggesting that both are real. Lower left: JVLA CO /= 1-0 imaging superimposed on our HST WFC3 Fl lOW snapshot, wherein 
we see faint diffraction spikes from the stars. Distorted, low-surface-brightness emission is seen coincident with galaxy C, characteristic of an interaction. Lower 
right: ~ 1.25^^ x 1.25^^ stamps centered on M, C, W and T, showing our JVLA 5.8-cm continuum imaging as blue contours superimposed on the HST FllOW 
imaging. At the high spatial resolution available here via interferometric imaging, we see a consistent morphological picture for the emission from cool gas and 
dust in the HATLAS J084933 system. The emission is dominated by W and T, which are separated from one another by ~85kpc in the plane of the sky, with 
significant emission also visible from two fainter galaxies, C and M. Individually, W and T are well resolved spatially. In all panels, contours are plotted at 
-3, 3 X cr, with v^-spaced increments thereafter, where a is the local noise level. Beam FWHM are shown as solid ellipses. N is up; E to the left; offsets from 
ck;2000 = 132.3 8 89° , (52ooo = 2.2457° are marked in arcseconds. Known stars are labelled. 



6 



IVISON ET AL. 



Wavelength / m [observed] 
10"' 10"^ 10"^ 10"'' 10"^ 10"^ 




Frequency / Hz [observed] 



Figure 2. Radio-through-optical SEDs of the dominant components of HAT- 
LAS J084933, W (red points, solid line) and T (blue points, dotted line). Our 
measurements cover the SED peak, constraining Ljr and Tdust well. The 
SEDs can be fitted adequately using a combination of thermal dust and syn- 
chrotron emission c omp onents (see ^3.3) plus the stellar population mod- 
elling described in 33.41 The key resulting characteristics, Ljr, T^ust, ^dust^ 
Mstars and ^ir are listed in Table |2] 

Streams, but rather significant concentrations of stars, gas and 
dust - luminous SMGs in their own right. 

Our deep rest- frame 133-289-nm Keck spectroscopy, with 
sHts covering W, T and C, reveals no significant line emission 
from T and C; the only line visible from W, albeit very faint, 
isCiiat232.6nm. 

Our CO /=3-2 image from CARMA is shown in Fig.[I] It 
is interesting to note that the low-level southern extension to 
the J =3-2 emission from T is mirrored at 3.6-4.5-/im, which 
suggests that both are real. 

3.2. Lensing model 

Galaxy W suffers no significant lensing magnification, as 
far as it is possible to discern from our extensive panchromatic 
imaging. 

The lenticular galaxy that lies along the line of sight to 
component T lies at z = 0.3 478, as revealed by our optical 
spectroscopic observations ( ^2. St , where the slit was placed 
along the major axis of the lenticular. Correcting for the in- 
strumental resolution, the lines due to Mg and Ca H and K 
suggest the foreground galaxy has a velocity dispersion of 
190±80]cms-i. Its stellar mass is (1.1 ±0.4) x 10^^ M©, 
adopting the approach outlined later in ^3.41 We have mod- 
elled this galaxy with an elliptical = 0.3) mass profile, trun- 
cated smoothly at a radius of 25 Icpc (the exact choice has little 
effect since the lensed image is well within this radius). The 
line of sight to the z = 2.41 starburst lies close to the semi- 
major axis of this foreground structure and suffers a magnifi- 
cation of /i = 1.5 ± 0.2, where the uncertainty was calculated 
by varying the model parameters - orientation, ellipticity and 
velocity dispersion - by 10%. 

3.3. Spectral energy distributions of W and T 

Fig. [2] shows the SED of HATLAS J084933, concentrat- 
ing on the brightest components, W and T. Our SHARC-2 
350-/im imaging provides evidence that the relative SPIRE 
contributions of components W and T follow the average ra- 
tio seen in the resolved imaging, 1.22 ±0.07, as described 
in ^3.11 Perhaps unsurprisingly, then, we then find they 



V /km s ^ (for z = 2.4100) 
-2000 -1000 1000 2000 




-2000 -1000 1000 



V /km (for z = 2.4100) 

Figures. CO spectra of HATLAS J084933. Upper two panels: CO J =1-0 
spectra of galaxies W and T superimposed on solid representations of the CO 
J =4-3 spectra (the former scaled by 16 x to put them on the same 7], scale). 
The two galaxies both have extraordinarily broad CO lines but the profiles 
are noticably different, so they are not lensed images of the same galaxy. The 
profiles of the CO transitions are indistinguishable for each individual com- 
ponent; their 7b ratios are also similar (Table |2). Middle two panels: CO 
J =4-3 spectra of galaxies M and C, with line centers significantly redwards 
of W and T. Lower panel: JVLA CO /= 1-0 spectra of W and T along- 
side the significantly different profile from GBT/Zpectrometer (Harris etdH 
2012). Emission at velocities redward of W and T - from M and C and sev- 
eral other faint clumps seen in the JVLA data, spanning a ~100-kpc region - 
can account for the differences between the GBT and JVLA line profiles. 



contribute in roughly the same ratio, 1.55 ±0.44, to an im- 
mense overall bolometric luminosity of HATLAS J08493 3, 
LiR = (4.8 ± 1.0) X 10^^ L0 (Table [3, having corrected the 
contribution from T for gravitational amplification (§ l3.2t . 

The IR-through-radio photometry can be described ade- 
quately with a model comprising synchrotron and thermal 
dust emission, following Kovacs et al. ( 2010). Table [21 sum- 
marises a simultaneous fit to the dominant col d dust tem- 
peratu re, Tdust, to LiR and ^ir (as defined by iHelou et aH 
I1985L but where 5'ir is measured across Arest = 8-1000 /im, 
as is Lir), and to Mdust (for a characteristic p hoton cross- 
section to mass ratio, k^so = 0.077 m^ kg~^ - iDunne et"an 
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Table 2 

Properties of HATLAS J084933.4+021443 



W T M C 



08:49:33.59 


r\o Ac\ 

08:49:32.96 


r\Ci AC\ oo o/^ 

08:49:33.80 


r\o AC\ 00 r\ 1 

08:49:33.91 


+U2: 14:44. o 


1 AO . 1 /I .in n 
+1)2:14:39. / 


1 AO . 1 /I . /I C 

+(JZ:14:4j.o 


1 AO . 1 /I . /I C A 

+UZ: 14:4j.U 


Z.9 ± U.4KpC 


3.5 ± U.OKpC 


b 


b 


/.U ± Z.i KpC 


O.Z ± 1.4 KpC 


b 


b 


1 'J ^9 -u n ozi 

i J. Hz U.U'4 


i J>. iO ZlZ U.UJ 


19 Q _|_ A 9 
iZ.!7 ZIZ U.Z 


19 _|_ A 9 
iZ.O ZIZ u.z 




i JUU 






on o _|_ 1 A 

HZ i .O 


1 _|_ 1 1 

JO. i ZIZ 1 . i 


b 


b 


Z. zt U.Uo 


_|_ A 1 Q 
Z.J J zt U. 1 J 


b 


b 


n _l_ A AC 

9.32 ± U.U5 


A 1 A _l_ A ACa 

9.1U±U.U5'' 


b 


b 


0.49 ±0.06 


0.56 ±0.07 


A AC7 _L A A1 1 
U.Uj / ± U.U13 


A ATA 1 A A1 A 

U.U/9 ± U.U14 


138±17 


157 ±20 


16.0±3.6 




6.72 ±0.82 


7.65 ±0.96 


0.779 ±0.178 


0.108±0.019 


825 ±115 


610±55 


320 ± 70 


250 ± 100 


2.4066 ±0.0006 


2.4090 ±0.0003 


/I 1 n _L A AAA/I 

2.41 /o ± (J.U(JU4 


/I 1 1 _L A AAAI 

Z.4135 ± U.UUUi 


3cr < 0.24'^ 


3cr < 0.24^^ 


3cr < 0.24'^ 


3cr < 0.24^^ 


4.08 ±0.92 


5.73±1.13 


0.66 ±0.27 


1.16±0.37 


830 ±155 


555 ±75 


260 ± 95 


620 ± 160 


2.4077 ±0.0009 


2.4096 ±0.0004 


2.4173 ±0.0006 




3.76 ±0.27 


5.10±0.37 


0.66 ±0.10 


0.95 ±0.15 


985 ±115 


545 ±30 


320 ±50 


450 ±90 


2.4068 ±0.0002 


2.4090 ±0.0002 


2.4178 ±0.0003 


2.4149 ±0.0004 


10.5 ±3.9 


b 


b 


160 ±60 


b 


b 


0.93 ±0.24 


1.14±0.27 


1.29 ±0.60 


1.63 ±0.59 


0.48 ±0.07 


0.57 ±0.08 


0.72 ±0.20 


0.75±0.18 


350 ±10 


220 ±15 


20 ±10 


80 ±20 


7.8±1.0 


6.8 ±0.8^ 


b 


b 


11.51 ±0.05 


11.13±0.05^ 


11.11±0.15 


11.11±0.20 


11.04 ±0.05 


10.92 ±0.06^ 


10.11 ±0.09 


10.25 ±0.07 


370 


210^ 


760 


440 


11.38±0.12 


11.01 ±0.12^ 


;^10 


10.36±0.18 



Property 

R.A. (J2000) 
Dec. (J2000) 

Mean FWHM at 5.9 cm, 870 /xm and in CO /=4-3 
FWHM in CO / = 1-0 /kpc 
log LiR /L0 

SFR /M0 yr-i ( Kennicutt 1998 . Chabrier 2003 IMF) 

^dust /K 
^IR 

log Mdust /M© 
CO/=l-0/co/Jykms-i 
CO / = 1-0 L^Q /lO'^ Kkms-i pc^ 
CO /= 1-0 Leo /lO^L© 
CO/=l-0FWHM/kms-i 

CO/=1-OZLSR 

i^C^^O and ^^Q}^0 /= 1-0 Iqq /Jy km s-^ 
CO/=3-2/co /Jykms-i 
CO /=3-2 FWHM /kms-i 
CO/=3-2 ZLSR 
CO/=4-3/co/Jykms-i 
CO 7=4-3 FWHM /km s'^ 
CO 7=4-3 ZLSR 
[Cii]/[cii] /lO-i^Wm-2 
[Cii] L[cii] /lO^Lo 

^3-2/1-0 
^4-3/1-0 
Vmax /km 

7? /kpc 

log Mdyn /M0 

log MH2+He /M© , for OLQQ =0.8 M© (K km pc^)~^ 
SFE /M0 

log Mstars /M0 

^ Corrected for /i = 1.5 ± 0.2. 
^ Insufficient S/N to make a useful measurement. 
^ Adopting 5totai/Vak as seen for ^^C^^O. 
Fixed. 

|2000"; Du nne & Eales|[20Qll - where we fixed the frequency 
dependence of the dust emissivity, /3, to be 2.0, which is 
physically plausible and consistent with the best-fit value, 
2.08 ± 0. 15). We fixed the synchrotron power-law index, a = 
-0.75, where Sy(x.v^, and we also fixed 7 = 7.2 - the power- 
law index of a dust temperature distribution appropriate for 
local starbursts, <iMdust/<^^dust r^"^^, which is designed to 
offer a physically motivated treatment of the Wien side of the 
thermal emission spectrum. Both W and T have Tdust ~ 36- 
40 K, commensurate with similarly luminous, dusty starbursts 
at z ~ 2-3, when calculated using a power-law temperature 
distribution (MagnelH et al. 2012). 

We are sensitive to AGN activity via our optical imaging 
and spectroscopy and our mid-IR/radio imaging. We have 
found no evidence^ that an AGN contributes significantly to 
the LiR of HATLAS J084933, though of course we have no 
evidence of absence since even a powerful AGN can be hid- 
den very effectively by the quantities of dust-rich gas in this 
system. 

If contributions to Ljr from AGN are small, then W and T 
are generating stars at a rate of ^34 00 and ^ 1 500 M ^ t) yr~^ , 
respectivel y, for a Chabrieij (120031) IMF (iKennicuttI 119981 
though see Dwek et al. 201 1). The estimates of Lir for C and 
M given in Table [2] are scaled from their submm photometry 

Observations of higher-/ CO lines would be useful in terms of AGN 
diagnostics. 



relative to W and T, assuming the same SED shape. This sug- 
gests that both are close to being HyLIRGs in their own right, 
^7x 10^2 

Before closing our SED discussion, we add a word of cau- 
tion regarding our estimated IR luminosities. If W and T lie 
in the core of a proto-cluster, as we shall see in what follows, 
then lower-luminosity star-forming cluster galaxies may con- 
tribute to the flux densities measu red by the relatively large 
SPIRE beam at 250-500 /im (see NegreUo etaH 120051) . Al- 
though the calibration cannot compare with that of HerscheU 
SHARC-2 recovers less flux than SPIRE at 350 /im, a possi- 
ble manifestation of the proto-cluster. 

3.4. Stellar masses 

The rest-frame ultraviolet-through-near-IR SEDs of W and 
T are illustrated in Fig.|2l alongside the far-IR-submm-radio 
data. Flux densities are listed in in Table [T] Our high-quality 
HST, VIKING and IRAC images at 0.88-4.5 /im are compro- 
mised to a certain extent by the unfortunate super-position of a 
star and a lenticular galaxy near W and T, respectively, but the 
considerable uncertainties associated with estimating stellar 
mass dominate over the photometric uncertai nties (see, e.g., 
iHainline et al.ll2011l:[Michalowski et al.ll2012l) . 

For components W, T and C we find best-fitting monochro- 
matic rest- frame //-band luminosities, logLn, of 1 1 .82, 11 .46 
and 11.12L0, respectively. The implied stellar masses 
are logMstars = 11-38, 11.01 and 10.36 M© for W, T and 
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C, using MAGPHYS with the most recent vers i on of the 
stellar populati on models of iBruzual & CharlotI (|2003|) . a 
IChabrieri d2003l) IMF and the SEP prio r distr ibutions de- 
scribed by Ida Cunha. Chariot. & Elba3 (l2008h . following 
iRowlands et al. (2012). 

Emission from M can be seen at 3.6-4.5 fim but the photo- 
metric uncertainties due to blending with C and W are large 
so we state only that its stellar mass must be considerable, 
^ lO^^M©. 

There is little sign of the power-law SEPs that can betray 
the presence of AG N, su pporting the lack of AGN spectral 
features reported in §13.11 though of course this does not rule 
out the presence of a deeply buried, highly obscured AGN. 

3.5. Spectral-line characteristics 

Fig. p shows the JVLA CO 1-0 spectra of W and T 
alongside our IRAM CO /=4-3 spectra. These spectral-line 
data bring us velocity information for the first time, enabling 
us to see that although W and T do lie at approximately the 
same redshift and share similar line profiles and 7b ratios, they 
have significantly different line widths (Table O. Each one is 
thus a distinct, extraordinarily luminous starb urst. While the 
CO li ne profile of T is typical of SMGs (e.g. ^Bothwell et al.l 
I2013h . that of W is broader than most, with ?^ 1000 km s"^ 
FWZL Both profiles are flat-topped or perhaps double-peaked; 
they can not be described well by a single Gaussian. 

The CO / = 1-0 fluxes of W and T measured in our JVLA 
data sum to /coi-o = 1.05 ± 0.09 Jykm s"^ At first sight 
this is consistent with the iHarris et aH GBT measurement, 
1.04 ± 0.37 Jy km s"^. However, the CO /=l-0 line profiles 
observed by JVLA and GBT appear noticeably different in 
Fig. [3) emission in the GBT spectrum continues significantly 
redward (by ^500kms~^) of the emission that can attributed 
to W and T. 

Looking at our IRAM PdBI imaging - see Fig. [5]- after sub- 
tracting the significant continuum emission (several submm 
continuum flux densities are listed in Table [TJ, we see that 
both C and M are detected securely in CO /=4-3 (we also see 
3.6-4.5-/im, submm and radio continuum emission - recall 
Fig.[T]). These galaxies are therefore unquestionably real and 
they share approximately the same redshift as components W 
and T. Their positions are listed in Table |2l along with those of 
the two brighter components. Their CO emission is centered 
^500 km s~^ redward of W and T. 

As well as C and M, a number of other faint emission fea- 
tures can seen in the integrated CO /= 1-0 image (Fig. [T]). 
Summing C, M and these faint clumps leads to tolerable 
agreement between the GBT and JVLA line profiles. Al- 
though we must add 0.87 ± 0.19 Jy km s~^ to /^o^^' we must 
also increase /^q^ by ^15% to account for the ^8-per-cen t 
lower flux density assumed for 3C 286 by lHarris et al.l (120121) . 
and for losses due to attenuation of the faint clumps by the 
GBT primary beam. The agreemenQ between the total GBT 
and JVLA line intensities remains consistent to within ^ la 
and the profile measured by the FWHM ^ 22^' GBT primary 
beam is thereby reconciled with the JVLA data. 

3.6. Resolved imaging of the gas and dust 

Our interferometric images from JVLA, IRAM PdBI and 
the SMA have a spatial resolution of ~ 0.5'^ FWHM, or 

Note that the low-level GBT spectrometer baseline (at negative veloci- 
ties relative to HATLAS J084933) may have caused an underestimate of /^o^- 



~ 4kpc. They allow us to compare the sizes of the regions 
responsible for emission from the total molecular gas reser- 
voir (via CO /= 1-0 emission), from the star-forming gas (via 
CO /=4-3), from cool dust produced by supernovae (SNe) in 
regions of recent star formation (via submm continuum emis- 
sion) and from relativistic electrons spiralling in or near recent 
SNe remnants (via radio continuum emission). 

We find that emission from W and T is resolved in all these 
wavebands, with total flux densities significantly higher than 
the peak flux densities. In Table [2] we quote their deconvolved 
FWHlvH sizes at 5.9 cm, 870 jim and in CO / = 4-3 and / = 
1-0, as determined with 2-D Gaussian fits. 

In the absence of significant AGN-powered activity, ra- 
dio continuum observations trace regions where massive stars 
have recently been formed, with no obscuration. Of the trac- 
ers available to us, we might anticipate the radio emission 
to be the most compact, as indeed is measured. At 870 fim 
we see no evidence of a more extended emission component, 
heated by older stars - the deconvolved sizes of W and T are 
consistent, as measured in 5.9-cm and 870-/im continuum, as 
one might expect given the far-IR/radio correlation (Table [J]). 
The CO /=4-3 FWHM are slightly larger, but are consistent 
with the radio and submm values to within the uncertainties. 

The error- weighted average of the 5.9-cm and 870- /im con- 
tinuum and CO / = 4-3 FWHM measurements is (3.4 ± 
0.4)kpc X (2.4 ± 0.4) kpc (PA, 9 ± 15°) for W; for T, before 
correcting for the amplification, we measure (5.6 ± 0.7) kpc x 
(3.5±0.7)kpc (R\, 48± 16°). In CO /=l-0, however, W 
and T have considerably larger deconvolved sizes than those 
measured in th e submm or radio continuum, as seen for SMGs 
generally fe.g. llvison et al.|[2QTTI:lRiechers et al.ll2011|) . 

Summarising, we have powerful, direct evidence of star- 
bursts covering ^ 3-4 kpc FWHM, with still larger reservoirs 
of gas available to fuel future star formation, on scales of ^ 6- 
7 kpc FWHM. 

The large sizes measured in CO have been mentioned 
alongside the relatively high LrcTn/^iR ratios m easured to- 
wards a number of SMGs fe.g. lStacey et al.ll2bl0|) . Both have 
been taken as evidence that SMGs form stars across larger 
spatial scales than the compact, nuclear events seen in local 
ULIRGs. 

Our SPIRE FTS spectrum does not have sufficient spatial 
or spectral resolution to separate the sub-components, but we 
do detect the Pi/i fine-structure line of singly-ionised 

carbon, [C ii] 157.74 /im, albeit at low significance (2.7cr; Ta- 
ble [2]). There is also a ~2-cr feature close to the wavelength 
expected for [Ol] 63.18 /im, which could conceivably be as- 
sociated with the blue wing of galaxy W, but we mention this 
no further here. 

The well-known plot of L[cii]/Lir versus Lcq/^ir is shown 
in Fig. m having corrected Leo and L rcm by 2 x and 0.7 x, 
respectively, following Hailey-Dunshe ath et al.l ^20 10), to ac- 
count for the fact that the optically thick ^^C^^ Q /= 1-0 emis- 
sion i s visible only from the side facing us (iKaufman et al.l 
Il999h and that only ^70% of the [C ll] is thought to arise 
in photodissociation regions (PDRs). This plot is commonly 
used as a diagnostic of density, n, and FUV field strength. 
Go. We find values of ^ ^ 5 x 10^ cm"^ and Go ^ 800 for 
HATLAS J08493 3. This is consistent with the values seen in 

We use half-light radii to calculate the area appropriate for gas and star- 
formation surface densities, alongside a 0.5 x correction to the SFR or gas 
mass. We use FWHM simply as a convenient way to compare sizes in several 
wavebands. 
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Figure 4. Lrcm/Z^iR versus Lcq/Ar for HATLAS J084933, adapted 
from Stacev et al. ( 2010), with the Solar-metalHcity PDR calculations by 
iKaufman et al. (1999) for gas density, n, and FUV field strength, Gq, the 
physical parameters to which these two ratios are sensitive. Regions typ- 
ically occupied by Galactic star-formation regions, local ultraluminous IR 
galaxies (UL IRGs), normal galaxies and starburst nuclei are illustrated (e.g. 
IStacev et al ' 1991; Malhotra et al. 2001; Luhman et al. 2003); the dashed line 
shows the i-[cii]/^co ~ 4100 typical for starbursts and Galactic star-forming 
regions. HATLAS J084933 occupies the region consistent with starburst nu- 
clei. 

nearby starburst nuclei, as has been the case for the majority 
of the other high-redshift SMGs detected in [C ll], supporting 
the idea that they are scaled-up versions of local starbursts. 

3.7. Dynamical characteristics 

The most remarkable characteristics of HATLAS J08 4933 
are revealed by the superlative spatial resolution (§ I3.6K ve- 
locity resolution and sensitivity of our JVLA CO / = 1-0 and 
IRAM PdBI CO /=4-3 data, as shown in Fig. [51 

To measure the 2-D velocity structure, we fit the CO emis- 
sion lines, pixel by pixel, in each of the cubes. Initially, we 
attempt to identify an emission line in 0.2'' x 0.2'' regions, in- 
creasing this slowly until we reach S/N > 5, at which point we 
fit the CO profile, allowing the intensity, centroid and width 
to vary. Uncertainties are calculated by perturbing each pa- 
rameter independently, allowing the remaining parameters to 
find their optimum values, until Ax^ = 1 is reached. 

We model the two-dimensional velocity fields of W and 
T by constructing and fitting two-dimensional disk mod- 
els using an arctan function, v(r) = 27r~^ Vasy marc tan (r/r^), 
to describe the kinematics (where Vasym is the asymptotic 
inclination-corrected rotational velocity and r^ is th e effective 
radius at which the rotation curve turns over - see ICourteaul 
[1997). We fit six free parameters to the disk velocity field: 
Vasym, ^t, J ccntroid, PA and inclination to the sky. The x^y 
centroid and r^ must be within the field of view and the max- 
imum circular velocity must be < 2 x the maximum velocity 
seen in the data. Quoted uncertainties reflect the range of ac- 
ceptable models from all of the attempted model fits. 

Our velocity cubes reveal that the wide CO line that led 
to this system being selected as a candidate HyLIRG is due 
partly to the multitude of constituent galaxies and partly 
to large regular, systematic rotational velocities in the two 
brightest constituent galaxies. Both W and T are resolved 
into counter-rotating gas disks - the scenario that many stud- 
ies have predicted should lead to the m ost intense starbursts 
(e.g. .Mihos & Hernquist 1994 119961: Itaniguchi & Shioval 



fT998l; iBomeetiD 12003 lBemll2QQTl: iDi Matteo etan[2QQ7l: 
Salome et al. 2013). 

As shown in Fig.[6l we find acceptable fits for W and T with 
disk inclinations of 56 ± 10° and 49 ± 10°, at PAs of 26 and 
145°, with Vasym = 415 and 370kms~^ respectively. Compo- 
nent T is sheared slightly by the foreground lens, as evident 
from the velocity and a maps. Both galaxies show small-scale 
(30-60-kms~^ , r.m.s.) deviations from the best-fit model. The 
dynamical masse£3 of W and T are listed in Table [21 several 
X 10^^ M0, where Mdyn = Rv]^^J^\x?{i). Their velocity dis- 
persions, cr, corrected for the instrumental velocity resolution, 
are ^ 65 km s~^ . 

M has no strong velocity gradient. Its Vmax is 20kms~^ 
(40kms-\ peak to peak), a = 140 ± 15 km s'^ for both CO 
1-0 and /=4-3, and Mdyn - 1.3 X lO^^M©. ForC, in CO 
/ = 1-0 we find Vmax = 80kms~^ (160 km s~\ peak to peak), 
a=130±20kms-^ and Mdyn - 1.3 x 10^^ (in CO/=4-3 
we find slightly higher values for a and Mdyn). 

3.8. Star formation in disks off the main sequence 

For main- sequence galaxies, star formation is not confined 
to a compact starbursting nucleus (e.g. Elbaz et al. 2011); at 
z ~ 2.4, the specific SFR is approximately 2.8Gyr~^. W 
and T have sSFRs of 14.2 and 14.7 Gyr~\ so their starbursti- 
nessE3 sSFR/sSFRms ^ 5, placing th em w ell above the main 
sequence. As we saw in ^3.61 and ^3.71 - their star forma- 
tion is not centrally concentrated as would be the case in lo- 
cal starbursts. Instead, W and T are undergoing widespread 
star formation in disks that are supported partly by rota- 
tion, as envisaged for main-sequence galaxies, and share 
many of the properties of massive star-forming disks (e.g. 
Forster Schreiber et al. 2009). Some other factor - presum- 
ably their interaction - has taken them into the realm of sSFR 
reserved for starbursts. 

For the Q p a ramet eJEl (iToomrd 119641: 
iGoldreich & Lynden-Belj Il965|) . which describes the 
stability of disks and which has been found to tend towards 
unity in stable situati ons (in models with and without feed- 
back - Hopkins et all 120121), we find values of ~ 0.35 for 
both W and T, with a plausible range that falls just short of 
unity. This suggests the disks are unstable and that the gas 
therein is prone to condense and form stars on a timescale 
shorter than the rotational period. 

3.9. Gas properties and aco 

The brightness temperature {T\y) ratios measured for both 
W and T are higher than expected (Table [21), where we usu- 
ally see r^-7/1-0 - 0.52 ± 0.09 and r^_yy^ ^ 0.41 ±0.07 for 
SMG s (lH^iseLal.1 120101; ll vison et al.[|2QlTI: i BothweU et all 
|2013|) . which suggests that the thermalised CO J = 3-2 fill- 
ing factors may be close to unity and that we may have found 
similar sizes in CO /= 1-0 and/=3-2 if our spatial resolution 
in the latter had been 2 x higher. 

For a CO to H2+He conversion factor, aco = 0.8 
(Kkm s~^ pc^)~\ as commonly adopted for IR-luminous star- 
bursts where the g as is not in virialized individual clouds 
(iBolatto et al.ll2013|) . the gas masses we determine for W, T, 

R was taken to be twice the deconvolved half-light radius in CO / = 1-0. 

We cannot determine the //?8 starburst indicator (Elbaz et al .1 120 111) for 
HATLAS J084933 since we lack sufficiently deep rest-frame 8-/im imaging. 

Us ually ascribed to a seminal paper on the stability of stellar disks, 
IToomrd ([1964), which contains no mention of Q. 
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Figures. CO / = 1-0 and J = 4-3 imaging of HATLAS J084933. Left: Velocity-integrated continuum- subtracted CO emission from our JVLA and PdBI 
observations. In both panels, we mark the positions of the brightest four galaxies, C, M, W and T. Right: Dynamical maps of the galaxies in CO / = 1-0 and 
J = 4-3. To construct these velocity fields, we fit the CO emission in each cube at each pixel, recording the intensity, velocity and line width where S/N > 5. 
North is up; East is left. Inset: One-dimensional rotation curves for W and T, extracted along the major kinematic axis identified in the modeling. Both W and T 
have rotation curves which resemble disks; indeed, W and T are best described as counter-rotating disks, which iMihos & HemquistI (11996) predicted would lead 
to the most intense starbursts. 
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Figure 6. CO J = 4-3 imaging of galaxies W and T, showing left to right their intensities, resolved velocity profiles, line-of-sight velocity dispersions (cr), our 
best-fit kinematic model (see text) and the residuals. The color scale for the velocity fields is shown in each panel. We note that there is evidence for lensing 
shear in the velocity field and line-of-sight velocity dispersion for T. 
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C and M are 110, 83.7, 12.8 and 17.8 x 10^ M© (corrected for 
lensing). 

If we were instead to adopt the all- sample average 
gas/dust mass ra tio for nearby galaxies (Mgas/Mdust = 9lt.f^- 
ISandstrom et al.l20 13). we would determine gas masses for W 
and T of 19.0 and 1 1 .5 x 10^ M©, yielding gas/dust-motivated 
aco values of 1.4 and 1.1 (Kkms~^ pc^)~^ respectively. 
Due to the uncertainty in a^sso, etc., we regard these as ac- 
curate to a factor 2 x at best. For an environment with sub- 
Solar metallicity, for example, aco would rise in line with 

Mgas/Mdust. 

A dynamically-motivated estimate of aco can be made by 
determining the difference between the dynamical and stellar 
masses in each system, since this represents the total plausi- 
ble mass of gao- Via this method, we find aco = 0-6 and 
0.3 M0 (K km s~^ pc^)~^ for W and T, respectively, where the 
main source of uncertainty is the notoriously difficult estimate 
of stellar mass, made even more difficult than usual here by 
the possibility of photometric contamination by nearby ob- 
jects and/or buried AGN(s). 

Finally, we can estimate aco using its relationship with 
gas-phas e metallicity and CO surface brightness, as devel- 
oped by [Narayanan et al.l ([2012b ). exploiting the resolved 
imaging in CO / = 1- acquired here. For Solar metallicity 
dChapman et al.""2005) and the size s tha t have been measured 
accurately in our JVLA imaging ( ^3.6[ Table O we find aco 
= 0.8 and 0.7 M© (K km s"^ pc^)"^ for W and T, respectively. 

Our mean is aco = 0.8 ± 0.4 (Kkms~^ pc^)~\ in line 
with the value usually adopted for starburst galaxies, with a 
range similar to that reported by iDownes & Solomo n* ('1998'), 
0.3-1.3 M0 (Kkms~^ pc^)~^ To reconcile our dynamical 
measuremen t of ac o with 3.2 (Kkms~^ pc^)~\ which 
iGenzel et"an (120101) argue is appropriate for all disks, regard- 
less of redshift, would require a buried AGN to dominate the 
near-IR light, causing us to over-estimate Mstars- 

Adopting our best estimate for aco, we find that W and 
T are both approaching the accepted criterion for maximal 
starbursts, with star-formation efficiencies, LjR/Mgas, of 370 
and 210 L0 Mq . C and M are higher still, though with con- 
siderable uncertainties. The Eddington limit is believed to 
be around 500 Mq ; above this level, the radiation pres- 
sure should quickly expel the gas via radia tion pressure from 
young, massive stars on dust (Scoville"2004). 

Feedback due to the aforementioned radiation pressure 
will lengthen the total duration of star-formation episodes, 
or lead to a series of short bursts. Ignoring this, the gas- 
consumption timescale implied for the galaxies comprising 
HATLAS J084933 is Tgas M^JSFR ^ 15-45 Myr, consis- 
tent \ yith an interaction-driven period of en hanced star forma- 
tion (^ Scoville. Sanders. & Clemenslll986l) . as already hinted 
by the morphologically tortured state of C 

The stellar masses determined in ^3.41 for W, T and 

C suggest gas fractions, /gas = Mgas/(Mgas + Mstars = 36- 

45%. These gas fractions are the top end of the range that 
iNaravanan. B othwell. & Davd (|2012a) showed to be difficult 
to reconcile with cosmological simulations. 

3.10. Relevance for star-formation laws 

Data of the quality gathered here have the potential to test 
whether the rules and prescriptions developed to describe star 

Bearing in mind that any AGN contribution to the rest-frame near-IR 
luminosity will have led to an over-estimate of stellar mass. 



formation in the local Universe (^e. g. ISchm idtll 1 9591: iKennicutl 
|1989|) can be applied to to the interstellar medium of distant, 
gas-rich starbursts like HATLAS J08493 3. Are the data con- 
sistent with a simple, volumetric star-formation law in which 
the SFR is ^1% of the molecula r gas mass per local free - 
fall time, as argued by Krumhol z. Dekel. & McKe^ (l2012h ? 
Are their SFRs more closely related to the orbital periods of 
their e ntire galactic disks (^e.g. lGenzel et al.ll2010l:lDaddi et al.l 
1201 Q[). or to the quantity of gas above some d ensity threshold 
(^e.g.'Lada e t al.ll2010l:lHdderman et al.l l2010')? 

We begin this discussion in the upper panel of Fig. Ui a 
luminosity-luminosity plot of the observables, Ljr and L^Q' 
for low- and high-redshift star-forming galaxies (SFGs), e.g. 
z = disks, LIRGs and ULIRG s and z = 1-2.5 BzK galaxies, 
taken from' Genzel et al.l (120101) . and for SMGs. For this plot 
we need make no assumption about the appropriate value of 
aco- The quality and depth of our continuum and CO data 
mean that the errors bars for W and T can barely be discerned. 
Using L'f^Q measurements made only in CO /= 1-0, thereby 
avoiding the need to adopt a 7b ratio, alongside self-consistent 
measurements of Ljr, we find no compelling evidence that 
different relations apply to SMGs and other SFGs, in line with 
the findings of Ivison et al. ( 201 1). Here, we use the sample of 
SFGs from lGenzel et al.l (120 10|) and find essentially the same 
result. 

It is noteworthy that while component T contributes most to 
Ico or L'f^Q in all three observed CO transitions, component W 
is significantly brighter in the continuum at IR-through-radio 
wavelengths. We find a similar variation between W and T 
in relation to the far-IR/radio flux ratio, with ^ir = 2.30 ± 
0.08 and 2.53 ±0.13, respectively, spanning (and consistent 
wit h) the tight correl ation for star-forming galaxies (^ir ^ 2.4, 
e.g. lYun et al.ll200l]) . One might think their different line-to- 
continuum ratios would lead to a relatively large difference 
in their positions on the plot of Ljr versus L^q (Fig.|7l upper 
panel). In fact, W and T both lie comfortably amongst the 
scatter of SMGs, suggesting that this scatter may be intrinsic 
rather than due to measurement error. 

Moving to the standard surface-density Schmidt- Kennicutt 
(S-K) relation, X^sfr ^ Xlgas' shown as the middle panel of 
Fig. U\ it is now necessary to adopt value(s) for q^cq appro- 
priate for SFGs and SMGs. Genzel et al.l (12010 1) argue that 
3.2 M0 (K km s~^ pc^)~^ is appropriate for disks, regardless of 
redshift; our dynamical aco constraints preclude the adoption 
of this value for W and T. We use 3.2 and 0.8 M© (Kkms"^ 
pc^)~^ for SFGs and SMGs, res pectively, acknowledgin g that 
this is overly simplistic^ (e.g. Narayanan et al."2012b'). We 
find W and T nestled amongst the other SMGs again (Fig.|7]), 
but this time the SMGs are offseQ significantly from the re- 
lation for other SFGs. 

The lower panel of Fig. |7] shows the Elmegreen-Silk (E- 
S) relation between star-formation surface density and the ra- 
tio of gas surface density and dynamical timescaleE3. Taken 
at face value, W and T exhibit extreme star-formation ef- 
ficiencies. Several other SMGs are equally extreme, con- 

Changing aco by a factor two (in opposite directions) for the SFGs and 
SMGs would bring them in line with one another. Because of the potential 
for a large systematic error in measurements of stellar mass , amongst other 
things, we view such a change as plausible. 

Any change in the area used to calculate the gas or SFR surface density 
results in a point shifting parallel to the power-law relation seen in this plot. 

We use the rotational timescale, calculated using the half-light radius 
and the velocity at that radius 
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Figure 7. Plots relating the rate or surface density of star formation with the 
available gas. Components W and T are shown with error bars in each plot. 
Circled red symbols denote high-redshift measurements made directly in CO 
/=l-0. Low- and high-redshift star-forming galaxies (SFGs), e.g. z = disks, 
LIRGs and ULIRGs and z = 1-2.5 BzK galaxies, are taken from Genzel et al. 
(^2010) for all three plots. Top: Ljr versus L^p, where we find no compelling 
evidence that SMGs follow a different relation than the other star-forming 
galaxies. Middle: Surface-density S-K relation, where W and T lie amongst 
hyperluminous SMGs such as SMM J02399-0136 and GN20, offset from 
the relation for other star-forming galaxies. Lower: Elmegreen-Silk relation 
between star-formation surface density and the ratio of gas surface density 
and dynamical timescale. W and T exhibit extreme star-formation efficien- 
cies, even relative to SMGs. The fraction of their available gas converted into 
stars per dynamical timescale, Cr > 0.1. 



verting > 10% of their available gas into stars per dynami- 
cal timescale (cr > 0.1). iKmmholz et al.l (1201 2h have argued 
in favour of a simple volumetric star-formation law, linking 
Galactic clouds, SMGs, and all scales between the two, with 
projection effects leading to the observed scatter. Here we 
have discovered two disks exhibiting star formation of such 
ferocity that they are capable of converting > 10% of their 
molecular gas into stars with each rotation. We require a num- 
ber of parameters to conspire together if we are to escape the 
conclusion that the gas in these disks does not obey the si mple 
volum etric star-formation law proposed by Krumholz et al.l 
(I20r2h . 

Before leaving the S-K and E-S laws, we note that al- 
though LiR is an effective trac er of Y^^pp in Eddington- l imited 
star-forming disks, Ballan tyne. Armour. & Indergaardl(l2013h 
have argued that velocity-integrated CO line intensity is a 
poor proxy for X^g^g. Resolved observations of the mid-IR ro- 
tational lines of H2 may provide the ultimate probe of the star- 
formation relations in high-redshi ft disks, requiring a far-I R 
space interferometer such as FIRI dHelmich & Ivisonll2008h . 

3.11. Future evolution of the system 

What does the HATLAS J084933 system look like in the 
present day? 

The extreme starbursts that brought this galaxy to our at- 
tention via Herschel are likely to have b een triggered by 
one or more interactions ("En gel et aP 120101) . Certainly, C is 
distorted, morphologically, and the distances between these 
galaxies are consistent with pe riods of intense star formation 
during merger simulations (e.g.lS pring el et al. 2005). 

However, without knowledge of the dark-matter halo they 
inhabit and the transverse velocities of W, T, C and M, it is 
impossible to be sure whether these galaxies are gravitation- 
ally bound, let alone whether they will merge. We find that 
the total energy, 1/2 + <l> (km s~^)^, of each of these systems 
is negative - indicating that the galaxies are bound to one an- 
other - if their dynamical masses represent as little as 50% of 
the total. This is confirmed by N-ho&y simulations, starting 
from the observed velocities and positions on the sky, and as- 
suming that the transverse velocities do not exceed those in 
the line of sight. If these galaxies inh abit an expected dark- 
matter halo of mass ^ 10^^ Mo (e.g. lAmblard"etan I20TTI: 
iHildebrandt et al.ll2013h . it is difficult to escape the conclu- 
sion that the system is bound. 

Taking another approach, w e searched the (500/ /z)^-Mpc^ 
Millenium Simulation at z = 2 (ISpringel et al.ll2005l) . In each 
friends-of-friends halo we identified all of the sub-halos, as 
well as the most massive sub-halo, which we take to be 
akin to W. We then computed the velocity and spatial off- 
set between the most massive sub-halo and all the other sub- 
halos and tracked the ^ 200 cases where there are three sub- 
halos within a radius of 300 kpc (co-moving) with velocities 
> 90, > 600 and > 800kms~^ in projection (i.e. taking into 
account likely line-of- sight distances and transverse veloci- 
ties). Those halos, typically ?^ 1.6 x 10^^ M© at z = 2, end 
up with a median mass of ^ 3.8;^^ 9 x 10^"^ M© in the present 
day (90% lower limit, 1.6 x IO^^Mq) - structures that are 
thought to have a space density of a few x 10~^ Mpc~^ (e.g. 
iReiprich & Bohr inger 2002t IVikhhnin et al. 2009). 

The scale over which we see the galaxies comprising HAT- 
LAS J084933, ^100 kpc, is similar to the core radius of rich 
clusters of galaxies, and the separation of W and T (in terms of 
kpc and km s~^) is typical of binary galaxies (e.g. lTumerl 19761: 
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Figure 8. L^q versus CO /= 1-0 FWHM line widths for SMGs with red- 
shifts, fluxes and magnifications from the literature (red squares - for de- 
tails, see Harris et al. 2012, from whence this figure has been adapted). The 
power-law fit to these data represents an estimate of the intrinsic line lumi- 
nosity versus line width relation for typical SMGs. Blue circles represent the 
lensed sample from Harris et al. (2012). The point for HATLAS J084933 has 
been replaced with those appropriate for components W and T. W has moved 
well within the scatter of the intrinsic relation, with a predicted amplification 
consistent with unity. Th e pre dicted amplification for T is consistent with the 
lens model presented in ^3.21 /i = 1.5 ± 0.2. 

ISchweize3ll98 7') and within a factor 2 x of the projected sep- 
aration of the two dominant galaxies in the Coma cluster0 
We conclude that such systems can survive for a long time 
and are not necessarily in the process of merging. The four 
starbursts in HATLAS J084933 may instead be signposting a 
proto-cluster, where the majority of members are likely well 
below the detection thresholds of our panchromatic imaging. 
It is ironic that deliberate attempts to select clusters at submm 
wavelengths (targeting radio galaxies and quasars, f or exam- 
ple - Stevens et al. 2003; D e Breuck etal . 2004; Greve et al.l 
120071: iPriddey et al...2008,) have met with only limited success, 
while serendip itous discoveries are becoming relatively com- 
monplace fe.g. lRiechers et al .1120101) . 

3.12. On the selection of HyLIRGs, and their rarity 

Having resolved HATLAS J084933 into its constituent 
parts, if we place W and T individually on the iHarris et aP 
plot of L^o versus FWHM line width (Fig.© we find that W 
has moved further within the scatter of the intrinsic relation, 
with a predicted amplification close to unity. The predicted 
amplification for T, /i 2, is consistent with the lens model 
presented in ^3.21 

We can estimate the space density of objects as rare as 
HATLAS J084933 by utilising the semi-empirical model of 
iHopkins et al.l (12010) and Havward et al.l (120131) . We start 
with halo mass functions and merger rates over the red- 
shift range of interest d erived from the Millenium Simula- 
tion (ISpringel et al.l l2005f) . assigning galaxies to each halo 
following the abundance- matching techniques described in 
IConroy & W echsler (2009) to yield the galaxy-galaxy merger 
rate as a function of mass, merger-mas s ratio and red- 
shift. Utilising the stellar mass functions of iMarchesini et al.l 
(120091). with extrap olations to higher redshifts following 
iFontana et al.l (120061) . we calculate the merger rate for all 

Coma (Abell 1656) is a B(inary)-type cluster, with two bright, super- 
giant galaxies, NGC 4874 and NGC 4889, close to each other in the cluster 
core (Rood & Sastrv 1971). 



galaxies with a stellar mass ratio equal to or greater than that 
of W and T, with one component at least as massive as W. 
Integrating the merger rates between z = 2-6 (t = 2.4 Gyr), 
where we are sensitive to such events via //-ATLAS, results 
in a space density of approximately 10~^Mpc~^, simil ar to 
the sp ace density p redict ed for extreme HyLIRGs (Lapi et alj 
1201 li recall also ^3.111) . For the typical gas-consumption 
timescale found in ^3.91 tsmg ~15-45Myr, we then expect 
to find ^3-9 such galaxies in the ^lOOdeg^ explored so far 
in //-ATLAS. The vast majority of bright Herschel galaxies 
have been shown to be lensed (Negrello et al. 2010), but the 
low end of this range is plausible. This suggests that the 
most luminous phase of these galaxies has a short duration, 
< 15Myr, and/or that such systems require some other rare 
property - counter-rotating disks being one possibility. Stel- 
lar mass functions at high redshift have significant uncertain- 
ties associated with them; nonetheless, this calculation em- 
phasises the rarity of major mergers such as that between W 
and T and illustrates the efficiency of the method employed 
here to find extreme HyLIRGs. The observed surface den- 
sity of these systems is such that even with short durations for 
their most luminous phases, they can end u p as a subset of the 
massive current-day clusters discussed in §13.111 

HATLAS J084933 was selected on the basis of its broad CO 
line profile as a low-magnification, intrinsically luminous sys- 
tem, so it is ironic that its wide line profile is in part due to the 
multiple nature of the source and that future searches for sin- 
gle maximal-luminosity galaxies may require the selection of 
less extreme profiles. However, fundamentally we conclude 
that the method explored here for the selection of the most 
extreme starburst events in the Universe shows promise. 

4. CONCLUSIONS 

1 . Exploiting the relationship between CO luminosity and 
line width determined for unlensed starbursts, we have 
identified and removed gravitationally lensed sources 
from the brightest galaxies found in the widest ex- 
tragalactic Herschel survey to yield a sample of in- 
trinsically luminous galaxies. Here, we report deep 
panchromatic follow-up observations of the best can- 
didate HyLIRG system, HATLAS J08493 3, with sub- 
arcsecond spectral and spatial resolution, which have 
led to the discovery of at least four starbursting galax- 
ies across a ?^100-kpc region at z = 2.41, each with a 
significant mass of stars. 

2. The two brightest galaxies, W and T, are separated by 
~85 kpc on the sky. W suffers no gravitational amplifi- 
cation. T is marginally lensed by a foreground lenticu- 
lar. 

3 . The panchromatic SEDs of W and T reveal that both are 
associated with high-luminosity events - HyLIRGs. If 
contributions to Ljr from AGN are small, and to date 
we have found no reason to suspect otherwise, then W 
and T have logLjR - 13.52 and 13.16, /^ust = 36-40 K, 
and are generating ?^3400 and ^ISOOM© yr~^ of stars, 
respectively. 

4. Two other gas-rich galaxies, M and C, detected 
nearby, make HATLAS J084933 reminiscent of the first 
SMG found and the best-studied HyLIRG to date, 
SMM J02399-0136, which comprises a complex galac- 
tic nursery with three massive stellar components, each 
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seen during a different evolutionary stage (Ferkinhoff 
et al., in preparation). The four galaxies known to com- 
prise HATLAS J084933 are most likely bound. A^-body 
simulations and comparison with the Millenium Simu- 
lation suggest that in the present day the system may 
resemble a B(inary)-type cluster of mass, ^ 10^^^ M©. 

5. We detect [Cii] from HATLAS J08493 3, albeit at low 
significance. Its flux relative to Leo and Ljr is indica- 
tive of n ^ 5 X lO'^cm"^ and Go ~ 800, values simi- 
lar to those inferred for nearby starburst nuclei, which 
supports the idea that the galaxies comprising HAT- 
LAS J084933 are scaled-up versions of local starbursts. 

6. Sub-arcsecond interferometric imaging reveals that the 
two brightest galaxies span ~3kpc FWHM in CO / = 
4-3, and in submm and radio continuum, and roughly 
double that in CO /= 1-0. Alongside the detection of 
[Cll], this supports a scenario involving a widespread 
burst of intense star formation. 

7. Exquisite 3-D imaging from JVLA and IRAM PdBI 
reveal counter-rotating gas disks in W and T, a sce- 
nario that has long been predicted to lead to the most 
intense bursts of star formation. This suggests that sim- 
ilarly luminous galaxies will often be found in pairs. 
Their modest velocity dispersions mean that the disks 
are prone to instabilities on many scales, even that of 
the disks themselves; they are undergoing coherent, ex- 
treme starbursts, at close to the Eddington limit, which 
are not confined to their nuclei. Despite their disk- 
like morphologies and dynamics, their specific star- 
formation rates place them ~ 5 x above the main se- 
quence. 

8. W and T have CO line intensities and widths typical 
of the brightest SMGs, with slightly higher 7b ratios. 
Three independent estimates of the CO-to-H2 conver- 
sion factor - exploiting the dynamical mass, the mass 
of dust and a recent metallicity-dependent relation - 
suggest aco = 0.8 ±0.4 M0 (Kkms~^ pc^)~^ and hence 
gas masses in W and T of ^ 10^^ Mq. Their gas frac- 
tions, ~ 40%, are difficult to reconcile with cosmolog- 
ical simulations. 

9. We have determined that the single-dish CO line profile 
of HATLAS J084933 is exceptionally broad because W 
an d T have rela tively wide lines, broadened further in 
the [Harris et al.l (^012 ) integrated GBT spectrum by the 
modest velocity offset between them, and by contribu- 
tions from a number of faint CO clumps slightly red- 
ward of the two dominant components. We can thus 
reconcile the GBT CO line profile with that of the 
various galaxies found to comprise HATLAS J084933. 
These situations likely happen only rarely, even in the 
widest extragalactic Herschel surveys. We estimate that 
only a few systems of this nature will be found in every 
100 deg^ surveyed. 

10. Placing W and T individually on the iHarris et al.l plot 
of L^o versus FWHM CO line width, W and T move 
further within the scatter of the intrinsic relation, with a 
predicted amplification close to unity. We conclude that 
the method outlined here for the selection of extreme 
starbursts can indeed separate HyLIRGs from the more 



numerous, less luminous, strongly lensed population of 
IR-bright galaxies, and may also pinpoint some distant 
clusters of starbursting proto-ellipticals. 

1 1 . Plotting LiR versus L^q , W and T lie amongst the SMGs 
on a single relation that is consistent with all star- 
forming galaxies. They also lie amongst the SMGs 
on the classic Schmidt- Kennicutt and Elmegreen-Silk 
plots, but the SMGs now follow a relation that is dis- 
tinct from those of other star-forming galaxies. Their 
star-formation efficiency is significantly higher, and is 
difficult to reconcile with a simple volumetric star- 
formation law in which ~1% of the available gas is 
converted into stars during each local free-fall time. 



Understanding systems as luminous as HATLAS J08493 3 
represents a challenge. However, the observations presented 
here have shed considerable light on the characteristics that 
likely lead to the formation - however briefly - of the most 
luminous star-forming galaxies in the Universe. Observa- 
tions scheduled for Cycle 1 with the Atacama Large Millime- 
ter/Submillimeter Array promise to reveal more about this 
dramatic and complex proto-galactic environment, probing 
the structure of the disks and sensitive to the presence and 
influence of buried AGN. 
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